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Abstract 

We investigate domain wall formation in late-time phase transitions. We 
find that as in the invisible- axion-domain- wall phenomenon, thermal effects 
alone are insufficient to drive different regions of the Universe to different 
parts of the disconnected vacuum manifold. This suggests that domain 
walls do not form unless either there is some supplemental (but perhaps 
not unreasonable) dynamics to localize the scalar field responsible for the 
phase transition to the low-temperature maximum (to an extraordinary 
precision) before the onset of the phase transition, or there is some non- 
thermal mechanism to produce large fluctuations in the scalar field. The 
fact that domain wall production is not a robust prediction of late-time 

transitions may suggest future directions in model building. 
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I. INTRODUCTION 


There has recently been interest in late-time phase transitions 1 as a possible solution 
to the large-scale structure problem. It has been proposed that the domain walls 
formed in a late-time (after decoupling) phase transition may seed the structure that we 
observe today, without conflicting with the well established isotropy of the microwave 
background radiation. (However, more detailed calculations show that this need not 
be the case. 2 ) 

Previously, it has been taken for granted that domain walls with interesting sizes 
(100 kpc and larger) can form in a late-time phase transition. In this work, we apply 
the results of a previous investigation into the dynamics of late-time phase transitions 3 
to examine the dependence upon initial conditions of the scalar field for the scenario of 
wall formation. We find that the formation of domain walls is not a robust prediction 
of models with late-time phase transitions. 

Before considering why domain walls might not form in a late- time phase transition, 
it is useful to recall why domain walls do form in standard phase transitions involving 
theories with disconnected vacuum manifolds. Consider a scenario for domain wall 
production in a simple model with a reflection symmetry ( <f> <-> —<}>), described by a 
theory with a single scalar field <f> and potential V(<t>) = (A/4 — <7 2 ) 2 . The zero- 
temperature vacuum manifold has degenerate vacua (<j>) = ±a. At temperatures above 
the critical temperature Tc = 2a, the ground state of the system is (</>) = 0. Fur- 
thermore, at temperatures much above the critical temperature, the (thermal) mass of 
the field, = AT 2 /8, increases with temperature. This has the effect of localizing 
the field about the high- temperature minimum. Below the critical temperature the 
high-temperature minimum becomes a local maximum, and since the field is already 
localized about (f> = 0, small thermal fluctuations will tend to push the field to one or 
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the other low-temperature vacuum states. The direction of the thermal kick will be 
random on scales set by the thermal correlation length. Thus, regions of the Universe 
will settle in different vacuum states, establishing a domain- wall network. 

We will refer to this picture of domain wall formation as the thermal mechanism, 
because it is thermal fluctuations that push the field to different vacua. We also note 
that thermal processes also play a role in positioning the field to a location where the 
thermal fluctuations can drive the field to different minima. 

We will see that a (possibly crucial) ingredient missing in the late-time transition 
model we study is a mechanism to localize the field in a position where thermal fluctua- 
tions will kick the field to different vacua. This problem arises as a result of the peculiar 
thermal properties of the field involved in the transition, essentially due to the fact that 
the weak coupling of the field results in very small thermal fluctuations. This leads 
to a scenario where the field will not evolve to a position where thermal fluctuations, 
which will be small, will drive the field to different low-temperature vacua. 

There are two ways around this problem. One could arrange for large fluctuations 
of super-horizon size in the scalar field well before the transition. In this case the 
wall network produced in the transition will reflect the initial field configuration, and 
there is no reason to expect the wall network to resemble the standard network as 
always assumed. A second possibility is to postulate that some dynamics, not part of 
the original model, acts to position the field correctly for the thermal fluctuations to 
operate. We will see that because the thermal fluctuations are so incredibly small, the 
positioning of the field must be done with great precision. 

Our discussions are based on the particle physics models for late-time phase transi- 
tion discussed in Ref. 3. The simplest model associated with phase transition consists 
of a single real scalar field with self interaction However, as pointed out in 

Ref. 3, in the context of a late-time phase transition, the more reasonable models to 


2 


consider axe the ones in which phase transition is driven not through <j> self interactions 
with the background, but rather by (f> interactions with a background of some other 
field t/>, perhaps a fermion. The essence is to have interactions weak enough to 
provide a late-time, soft-wall transition, but ip can have additional interactions that 
can establish the background by thermal interactions. 

This paper is divided into four parts. In Sec. II we calculate the thermal fluctuations 
of the scalar field involved in the transition. In Sec. Ill we consider the dynamics of 
the evolution of the scalar field in the early Universe. In the final section we conclude 
and discuss the prospects for wall production in late-time phase transitions. 

II. THERMAL FLUCTUATIONS IN SOFT- WALL MODELS 

A. YUKAWA MODEL 

First, let us consider the simple model of a single real scalar field <f> interacting with 
a fermion ip, with a potential: 

KoO, ip) = ^ (2.1) 

where the parameters mo, A 0 , and h are the unrenormalized mass and coupling con- 
stants. Following the analysis of Ref. 3, we assume that the fermion loops dominate, 
which is the case if h >■ A 0 and the boson tadpole is neglected. The renormalized, 
one-loop effective potential at finite temperature is 

W) = -jk J W) , ln^ + AV' T W. (2-2) 

where \x is an arbitrary mass scale related to the values of coupling constants, and 

A V t {4>) = -4^- dxx 2 In 1 + exp ( -\Jx 2 + m^/T 2 ) , (2.3) 
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with m = h<f>. 

At high temperature, AVr(<f> ) — +h 2 <f> 2 T 2 / 3. Hence the curvature at <j> = 0 (which is 
negative for T — 0) will be positive for high T, and <)> = 0 becomes the global minimum 
of the potential. Defining the critical temperature for the phase transition to be the 
temperature when d 2 V/d 2 <j> vanishes at <f> = 0, we have Tc ~ We want the 

phase transition to occur after recombination (but well before the present epoch) and 
to produce “soft” walls. These requirements lead to restrictions on ( m ^ < 10~ 24 eV) 
and Tc (eV > Tc > 10 -3 eV). By making h sufficiently small, it is possible to have 
Tc for a late-time, soft-wall phase transition. 

Thermal fluctuations in the field are found by examining the finite-temperature, 
two-point correlation function. The two-point correlation function is given by the 
Green’s function 

(mm) T = gw *) = / ^ *- <k *D T (k) t ( 2 . 4 ) 

where Dr(k) is the <j) propagator at finite temperature. In keeping with the assumption 
that <f) is very weakly coupled, we will assume that it is not in equilibrium and there is 
no thermal (f> background. Thus the thermal effects enter only through the interactions 
of <j> with the thermal background of tps. To compute Dr(k), consider the one-loop 
diagram in Fig. 1. Recall that the fermion propagator at finite temperature is 4 

S T (p) = i(i> - m^r 1 - 2 tt + m^)6(p 2 - mj,), (2.5) 

where / \p is the phase space density for Let us assume that the phase space density 
for is that of a thermal distribution (i.e. , a Fermi-Dirac distribution with temperature 
T). With this we have 

D T (k ) = i(k 2 - mj) -1 + i(k 2 - mj) _1 • loop ■ i(k 2 - mj) _1 , (2.6) 

where 
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loop = (-1) J Tr { ihS T (p ) ihS T (p - *)} • (2.7) 

We will ignore the temperature-independent contribution of the loop to Dj(k), as 
it is absorbed in the renormalization of the (f> propagator, and concentrate on the 
temperature-dependent contribution. Denoting the temperature-dependent paxt of the 
loop as loop?, we have 

loopT = K - k)+ m *l s (p 2 - o 

+4 xV*(£p)/*(E p - t )«(P 2 - mJ)<S (0> - kf - m\) 1 . (2.8) 


The leading term in this expression corresponds to taking the thermal part of one 
fermion propagator and the temperature-independent part for the other, i.e., the first 
term. The temperature-dependent two-point correlation function is then given by 

g t(i) = J i*pU(E„)6(p 2 - mj) 

x / d ' k c -ifa Kp - *) + f2 91 
(2 tt) 4 (A: 2 - m^) 2 [(fc — p) 2 — 

Now at finite temperature there are three length scales in the problem that appear 
in Gt{x): m^ 1 , m^ 1 , and T -1 . It is useful to define the following dimensionless length 
scales: 


= W/m/; P* = M/m/; p T = |x|/T z = p^/p T . 


( 2 . 10 ) 


After much manipulation, we may express Gt{x) in terms of sums and integrals in- 
volving modified Bessel functions: 


gkixI) * 


f Jfl( ;^ +(1 -^) ^ (yd -m+'4) 

* In 2 + (1 — 9) 2 Pt 
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Ki (z^p + (i - ey£) 


+ *n)f»e*'y n ‘ 

Jo sjn‘ + (l-e) 


- m Vr 


( 211 ) 


\/( 1 - oy + dml/ml j ’ 
where we have used m $ <C rrty, which will be true in the context of late-time phase 
transitions. Note that for jxj = 0, we have 


<?r(0) = 


h 2 


Ki(nz), 


( 2 . 12 ) 


(2tt) 3 ^ n 

which for T m^, gives Gt( 0) = (h 2 /96ir)(m^/m,p)T 2 , the familiar T 2 dependence of 
thermal fluctuations in 4> 2 . 

For |x| ^ 0, we axe able to obtain an upper limit for Gr(|x|): 


£r(|x|) < 


7 ih 2 m^T 
(27r) 4 


A r i(p 0 )ln (l + e 2 ) , 


(2.13) 


again assuming that m ^ m^. 

For the analysis of domain wall formation, we will be interested in correlations in 
4> on length scales comparable to the scale of the correlation length in the transition, 
l x l ~ m 4 , l i i- e -j P4> ~ 1- Since Tc ~ m^/h, and Tc < 0.1 eV ~ m ^ as discussed 
previously, it follows that that h ~ Hence 


Gt(P4> — 1) < rn^T Ki(p^). (2-14) 

Since p^ <C pr for late-time phase transitions, it is no surprise that the relevant 
dimensionless length scale. is p Note that ~ 1) is of order unity, while Ki(p<t> 

1) ~ ^7r/2p^ exp(— p^). It is clear that the correlations in <$> are rather insubstantial 
on scales relevant to domain wall formation since C Tc ~ m^. 

This simple model for late-time phase transition involves a very small, perhaps 
unnatural dimensionless coupling constant h. In the next section we study a more 
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realistic model in which there is physical motivation for the origin of small coupling 
constants. 


B. CHIRAL MODELS 

Let us now consider a theory with the Lagrangian proposed by Hill and Ross: 5 


C = ^d^4>d^(f> + u L i $v L + u R i flvR + (rnv L URe x4l l s + eu L u R + h.c.'j 

+m 4 cos + <5^ . (2.15) 

Note that the above theory has non- linearly realized chiral U(l) symmetry when the 
e-term and the m-term are not present. The m- term has been added to give <f> a mass of 
unspecified origin, this usually comes from some deeper (compared to the spontaneous 
breaking of the original chiral U(l) symmetry) symmetry breaking in the theory which 
breaks the continuous U(l) down to a discrete subgroup 2^. The e-term can be 
interpreted as the additional manifestation of this deeper symmetry breaking, similar 
to the cr-term in the nucleon-pion system. It breaks the chiral symmetry explicitly. 

Since a nonzero m term will always be induced by the presence of a nonzero e and 
m, let us write down the following Lagrangian which contains N Dirac neutrino species 
and is invariant under a Zn discrete symmetry: 

C = Uj i fluj + XI ( m + ee’' A//+2 ‘- 77r/ ' ,v ) UjL^jR + h.c. (2.16) 

1 j = o j=o 

First consider the N = 2 case. The effective low-energy theory (energy scales much 
less than /) has Yukawa couplings of the form: 


C Y uk = M + (<f>)^ 0 ip 0 + M_(^,^i 

M|(^) = m 2 + e 2 ± 2mecos(<f>/ f). (2-17) 


It has been shown in Ref. (3) that the finite temperature potential for <p in the Z 2 
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case is 


Vt(4>) 

mi 


1 (1 + T'f In (^' 


327 r 2 

+ 


, ('• 2 ^ A )/' r2lnr 2 + 


2 mi 


+ 


r 2 -l' 


32tt 2 


21 f 

cos — 


m_ y mi \ m_ y mi 


2tt 2 

where we have defined 


(; 

^2 J x 2 In [l -f exp (-^Jx 2 + M/(<£)/T 2 )] (2.18) 

v / j— 4“ > 


' m± = m ± e; r = m + /m_; = V r// (</> = 0). 


(2.19) 


The corresponding temperature-dependent mass is defined as m^(T) = V^'(^> = 0): 

2 meT 2 r°° , •> (x 2 + m?/T 2 )~ 1/2 


= m o + 


.? 


da: x^ 


7i- 2 / 2 J ' Jo 1 + exp(x 2 + m^/T 2 ) 1 / 2 

2 2meT ^ (— l) n+1 / m_\ / m+\ 

= m ° + 7vr£— s— l m - K '( n l r )- m+jr, ( n l r ) 

For T m ± , we get 

T 


m 0 (T) = m 0 + 


me 

^T 2 


ml In 


m+ 


T 

mi In 

m_ 


( 2 . 20 ) 


( 2 . 21 ) 


Obviously, at high T, we have m^(T) > 0, i.e., we have a phase transition with Tc 
satisfying m^(Tc) = 0, provided that m^ < 0. 

By expanding the mass functions M±(<f>) in Eq. (2.17), we find the lowest-order 
coupling of 4> to the fermions: 


A = h + 4>4>tp 0 xpo + h^<jxf> t/qVh 

, 1 me 

h± = T- 


„ • ( 2 . 22 ) 

m± P 

To compute the two-point function, we need the lowest-order correction to the <f> 
propagator, see Fig. 2. The correlation function is still given by Eqs. (2.4) and (2.6), 
with 
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loop = (- 1 ) J Tr (ih+s£(p) + ih-S T (p ) ) , 


where 


(2.23) 


Sf = i(]i ~ m±) 1 - 2 Trft(E p )(rf + m±)S(p 2 — m\.). 


(2.24) 


Let the temperature-dependent part of loop be loop 7 . We then have 


loopT = loopj + loop T 


loop ± = 


(2?r) 2 

Eq. (2.4) now gives 


k±<T t i^f) 


n=l 


(2.25) 


£r(M) = xo -KoOv) • loop T . 


2(2t r 


Hence we obtain 


Gr( ' X|) = 2&7p K °M 

00 ^_2jn+l 


X 




(2.26) 


(2.27) 


For T m±, we have 

Gt(W) -( 2 ^f(t) k ° m ln(r/ro) ’ (2 ' 28) 

Note that C?t(|x|) is singular for |x| — ► 0, which is natural since the effective Lagrangian 
of Eq. (2.15) is not renormalizable. 

Gt(|x|) is well behaved for |x| 7 ^ 0 and gives us the correlation function between 
two spatially separated points in the <j> field. Note that 


Kq{P4> > 1 ) 


7T 


4> 


exp (- P +), 


(2.29) 


9 


while Ko(p<t> ~ 1) < 1. Obviously, the correlations in <f> field are suppressed severely on 
the scale of |x| ~ m^ 1 , since e < m <C /. 

In soft-boson particle physics models associated with late-time phase transitions 
we have the typical scales ~ O.leV; / ~ 10 15 GeV; Tc ~ and m# ~ m^/ f. 
The domain walls are cosmologically interesting because the walls have scale m^ 1 ~ 
m V’ 1 (// m V')- The smallness of (m^//) is essential in making the wall scale large enough 
to be of interest to large-scale structure, however, it also leads to suppression of the 
correlation in the <f> field due to thermal effects. Stochastic noise is washed out in the 
large volume m^ 3 . 

III. THE EVOLUTION OF <f> 

In the last section we demonstrated that a generic feature of the late-time, soft- 
wall models is very small thermal noise in the system. This means if wall formation 
follows the standard scenario, there must be some mechanism to poise the field at a low- 
temperature maximum to an extraordinarily high precision, one part in (//m) 2 ~ 10 5 °. 
Is this reasonable? Well, since the low-temperature maximum is a high-temperature 
minimum, the field naturally will evolve toward the desired point and undergo damped 
oscillation about that point. In this section we show that, as one might expect, the 
oscillations will not be damped enough to set the field to the desired point to the 
necessary precision. 

Let us now turn to the evolution of the (f> field in the potential Vr(<j>) given by 
Eq. (2.18). In Fig. 3, Vr(<f>) is plotted for several different temperatures. The phase 
transition occurs when the potential turns over at T = Tc- Solving Eq. (2.20) numeri- 
cally with m^(T) = 0, we obtain Tc as a function of (—mlf 2 /m 4 _), given r = m + /m_. 
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The result is plotted in Fig. 4. Note that the high temperature behaviour of the Tc 
curves agrees with the formula obtained from Eq. (2.21): 3 


Tc = m + 


m + 


m l/(m2 —m?_) 


exp 


' 8ir 2 f 2 (—miy 


(3-1) 


m.j - \(ml-ml) 2 )' 

The issue at hand is if the “initial” value of 4 i», <f>(T Tc) is displaced from the high- 
temperature minimum, will it relax sufficiently close to the high-temperature minimum 
by Tc so that it can receive random thermal kicks to different regions of the low- 
temperature vacuum manifold. To address this issue we follow the equation of motion 
for the scalar field: 


• OVr 

<f> + 3 H 4> H — ;rr~ = 0 


d<f> 


3 30 


H 2 = 


9*T\ 


(3.2) 


where we have assumed a radiation-dominated Universe. 
Using the dimensionless quantities 


c-t 

9 = r 


T = 


ra_ 


V tW = — 5-; 

ml 


7 = ~ mlP 


mi 


(3.3) 


after some manipulation, the equation of motion becomes 


d?9 


90 


f Mpl \ 2 1 dvr 


(3.4) 


dT 2 8* 3 9* \ f ) T 6 99 
There is no analytical solution to this equation. A numerical solution is shown by the 
solid curve of Fig. 5. As expected, the <f> field oscillates about the high- temperature 
minimum (here <j> = 0), then after the phase transition it evolves to a low-temperature 
minimum (here (j> = — 7t/ 2). The reason the field evolved to <f> = —tt/2 rather than to the 
degenerate low-temperature minimum (j) = +7r/2 is simply because at the temperature 
of the phase transition the field was negative in its oscillations about <f> = 0. 

In order for walls to form we want some mechanism to send the field to <f> = — -rr/2 in 
some regions and <f> = +7t/ 2 in other regions. Now if we wish we can model the thermal 
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noise in the system by the well studied stochastic approach. 6 However for thermal 
fluctuations to push cf> to different regions, they must have an amplitude comparable 
to the oscillation amplitude. We saw from the previous section that the amplitude 
of thermal fluctuations are very small indeed. Is there any way for the amplitude of 


oscillations about <f> = 0 to be damped sufficiently to empower thermal fluctuations? 


Clearly for the parameters of Fig. 5 the answer is no. 

We can explore the sensitivity of the oscillation amplitude at Tc upon model pa- 
rameters by considering the adiabatic invariant of the oscillating system. 7 
For oscillations near 9 = 0 


dv T J_d¥L _ JL m 2 (T\9 

de mi d(f> mi 

and we have cP6/dT 2 = —A 2 9, with 


^ = / 90 / Mpl\ m^{T) 


8tt 3 <7* \ mi 


(3.5) 


(3.6) 


If A were constant, the system would be closed and would execute a strictly periodic 
motion with a constant energy E and a fixed period A -1 , i.e., 9 = Acos(AT). When A 
is variable, the system is not closed and its energy not conserved. If A changes slowly, 
the rate of change in E will be small. The dependence of E on A can be expressed as 
the adiabatic invariant I = E / A ~ AA 2 , when the adiabatic condition \\~ l d\/ dT\ <C A 
is satisfied. 

The <}> field first begins to oscillate at t ~ m^ 1 . Since T 2 = y^O/87 r 3 gr* ( MpL/ 2 t ), 


Tosc = T(t - m^ 1 ) = 


/l / 90 \ 1 Im^MpL I Mpl 

UVSvJ V - mi 0(1 


(3.7) 


It is straightforward to show that at T osc we can use the adiabatic approximation. Of 
course at T — Tc , the adiabatic condition breaks down, because m <t> {Tc') = 0 leads 
to A(Tc) = 0. If we are only interested in a lower limit to the amplitude of the 
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oscillation, we may set m $ equal to a constant, say its zero- temperature value. In this 
case A cc A" 1 / 2 ~ T z >\ and A(T C < T < T osc ) ~ A(T osc ) x (T/T osc ) 3 / 2 . 

Now we can compute the amplitude of the classical oscillations near 9 = </>/ f = 0 
for T close to Tc (assuming Tc < T osc ; if the inequality is not obeyed, there is no 
damping): 

MT ~ T c ) =* A(T = r o «) (|hY , (3.8) 

N •* OSC / 

and use Eq. (3.1) for T<? and Eq. (3.7) for T osc . 

An example of this approximation is shown in Fig. 5. The solid curve is a result of a 
numerical integration of the equation of motion, and the dashed curve is the adiabatic 
approximation for the oscillation amplitude. The parameters chosen for the figure 
are f/Mpi = 10~ 4 , / = 1, m + = 3m_, and <j)(T Tc)/f = 1. We find that the 
adiabatic approximation holds for the range of parameters of interest. For the entire 
parameter space the amplitude of oscillation is never reduced by more than a few orders 
of magnitude. 

The conclusions are clear: There are insufficient oscillations between T osc and Tc to 
damp the amplitude of the oscillations to the desired scale of [see Eq. (2.28)] 

= sfmmF* ~ i)) T = \A?r(M) = ^ (y) vCt, (3.9) 

which would correspond to A(Tc ) ~ me/ f 2 ~ 10 -50 . 

Since thermal effects can not pin <f> to a value where thermal effects can kick the (f> 
field to different minima, <j> has to be initially localized (to extreme accuracy) to <f> = 0 
via some non-thermal mechanism for domain walls to form via the thermal process. 

Finally, we mention that although we have only studied the Z 2 model, our conclu- 
sions also hold for Z;v > 2 models as well. In fact, for these models the high-temperature 
<f> mass is driven to zero, so there is even less of a force driving <j> to the appropriate 
position. 


13 


IV. CONCLUSIONS AND PROSPECTUS 


We have studied the thermal effects on the two-point <f> correlation functions for 
two models associated with late-time phase transition. The first model is the simple 
A <p 4 theory with Yukawa coupling between <p and some thermal fermion ip, the second 
model is more natural and involves a cosine potential for <p induced by its coupling to 
thermal neutrinos ipo and ipi. We have shown that in both cases the thermal effects 
are suppressed due to the smallness of m $ compared to m^. Since the essential feature 
of a late-time phase transition model is to have domain walls with cosmologically large 
scales caused by the smallness of m^, we conclude that domain wall formation does not 
follow the usual thermal scenario unless there is non-thermal mechanism for localizing 
<p to zero initially. 

Of course one might imagine that the formation of domain walls in late-time phase 
transitions does not follow the thermal scenaxio, but results simply because long before 
Tc there are large, non-thermal fluctuations of the field, and the classical evolution of 
the field takes the field to different regions of the vacuum manifold. In some sense, the 
walls are present as initial conditions at Tc- If one adopts this approach, a domain wall 
network will result with characteristics (wall sizes, distributions, velocities, etc.) that 
reflect the initial conditions, and need not resemble the properties of walls produced by 
the thermal process. One must be careful to arrange these fiuctutions before inflation, 
since inflation will erase them on scales of interest. 

This non-thermal scenario is far from impossible, since the starting Lagrangian is 
of course only a low-energy approximation to the complete theory. It may well turn 
out to be natural that the seemingly highly contrived conditions necessary for domain 
wall formation are a natural result of the ultimate theory. Until such models are 
constructed, it is impossible to make predictions sensitive to the domain wall network, 
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since the network depends upon initial conditions. 

Our work suggests that if late-time phase transitions play a role in structure forma- 
tion, it is more likely to be through slow-roll scenarios like those studied by Frieman, 
Hill, and Watkins 8 than through formation of domain walls. 
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FIGURE CAPTIONS 


Figure 1. One- loop diagram of the <f> propagator for the model with simple Yukawa 
couplings. 

Figure 2. One-loop diagram of <f> propagator for the Z 2 soft-boson model. 

Figure 3. The total potential Vr(<j>) for the Z^ model, plotted for several different 
temperatures. The phase transition occurs when the potential turns over at T — Tq. 

Figure 4. Numerical solution to m^(Tc) = 0. Tc is plotted as a function of 
(— mo/ 2 /ml.), for several different values of r = m+/m_. Read from top to bottom, 
the curves are for r = 1.01, 1.05, 1.5, 2, 3, and 5, respectively. 

Figure 5. Example of the evolution of the <f> field. The solid curve is a result of 
numerical integration of the equation of motion, and the dashed line is the adiabatic 
approximation for the oscillation amplitude. Parameters for the model are f /Mpl = 
10~ 4 , / = 1, m + = 3 m_, and <f>(T ^ Tc)/ f = 1. 
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